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The Hubbard model of bosons on two dimensional lattices with a lowest flat band is discussed. 
In these systems there is a critical density, where the ground state is known exactly and 
can be represented as a charge density wave. Above this critical filling, depending on the 
lattice structure and the interaction strength, the additional particles are either delocalised 
and condensate in the ground state, or they form pairs. Pairs occur at strong interactions, 
e.g., on the chequerboard lattice. The general mechanism behind this phenomenon is discussed. 


1 Introduction 

Flat band systems have been studied extensively in experiment and theory. They are a prototype for 
strongly correlated systems. Strongly correlated phases of matter emerge in such systems since the inter¬ 
action, even if it is small, dominates the behaviour. The interaction leads to different phases. Fermions in 
flat bands are known to show ferromagnetism [H El [3] and ferrimagnetism H] , for a review we refer to [5] . 
More recently, bosons in flat bands have been studied, esp. the question of Bose condensation and the 
emergence of topologically ordered phases such as lattice versions of fractional quantum Hall states (see 
e.g. [ZlEllg]) are important here. Hard-core bosons can be mapped to spin systems, which as well have 
been investigated, see e.g. |Tnl flT] . 

Promising experiments to realise such systems are ultra-cold atoms in optical lattices; as they provide 
a perfect implementation of the Hubbard model. Optical lattices might be viewed as quantum simulators 
and are realised by coun-ter-propagating laser beams forming a periodic microtrap for the atoms. These 
experiments enable the control of a large number of parameters as the potential depth or the lattice 
geometry itself, see e.g. nainiiiiiis]. 

For bosons in a flat band, different phases may occur. A especially interesting question is whether 
Bose condensation occurs. This question was, to our knowledge, first investigated by Huber et al. for 
bosons on a kagome lattice. The kagome lattice is the line graph of the hexagonal lattice |Tj and thus 
a prototype of a larger class of lattices. For these and similar lattices, the dimensionality is important. 
This was already shown by Huber et al. |6]. They compared the one-dimensional saw-tooth chain, the 
line graph of a decorated chain, with the kagome lattice. Both lattices have an ordered ground state at 
some critical filling pc'- a charge density wave (CDW). Adding further particles to the CDW destroys the 
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order in the one-dimensional case. Domain walls appear, which can freely move, and the system becomes 
a Luttinger liquid |6]. This result holds at weak interaction. In the hard-core limit, a two-body bound 
state is formed which can move through the lattice nsunKiB]. On the kagome lattice, in contrast, domain 
walls are less favourable than extended states. The CDW ground state on the kagome lattice is three-fold 
degenerate. Each of these states has interstitial sites which are unoccupied. The additional particles have 
a high mobility and probability to move on the interstitial sites. A small but finite density of additional 
particles form a condensate. The underlying CDW states remain intact for densities above but close to the 
critical density and break down at higher densities [6]. These results hold for weakly interacting bosons. 
For strongly interacting bosons, esp. for hard-core bosons, other states may be favourable. 

Several questions arise from these results. Since the kagome lattice is a prototype of a larger class of 
lattices, one may ask if the results are true on other line graphs of two dimensional bipartite lattices. 
Another important member of this class is the chequerboard lattice. For the chequerboard lattice, the 
CDW state at the critical density is two-fold degenerate and there are no interstitial sites [19]. The first 
question therefore arises, whether on the chequerboard lattice the physics of weakly interacting bosons is 
the same as on the kagome lattice. A second question is whether for strongly interacting bosons on one of 
these lattices bound states occur, as in the one-dimensional system. A third question is, what happens in 
three dimensions. 

In this paper, we concentrate on line graphs of two dimensional bipartite lattices, for several reasons. 
On this class of lattices, the flat band contains strictly localised states. We will briefly sketch them and 
their properties in the next section. For details we refer to |19[ I20j . Each finite face of the two dimensional 
bipartite lattices is surrounded by a cycle with an even number of sites. In |2n| it was shown that the 
states which are localised on these cycles are linearly independent and complete. They form a basis, but 
not an orthonormal basis, which makes it difficult to deal with them. For that reason, we follow the idea of 
Huber et al. [^ and introduce a Wannier basis. The Wannier states are as well localised, but not strictly 
on few lattice sites. They fall of algebraically. Similar to the approach by Huber et al. we use the Wannier 
states in Sect. to derive an effective Hamiltonian for weakly interacting bosons on the chequerboard 
lattice. It turns out that the essential physics for weak interactions is the same as on the kagome lattice. 

Let us mention that the algebraic decay of the Wannier states was as well used by Chalker et al. |21] in 
their investigation of disorder and localisation in two-dimensional flat-band systems. 

In Sect. I^we study the case of hard-core bosons. We use different kinds of variational states, partly 
by using exact diagonalisations of small separable sub-units of the lattice. On the chequerboard lattice, it 
turns out that pairs of bosons are formed. On the kagome lattice, the pair states exist as well, but have not 
the lowest energy. The ground state on small separable sub-units of the lattice shows a larger expectation 
value of the particle number on the interstitial sites. This indicates that on the kagome lattice the physical 
picture developed by Huber et al. remains true for strong interactions. The physics of hard-core bosons on 
the two lattices, kagome and chequerboard, is thus different. Finally, in Sect. |^we give some conclusions 
and discuss open questions. 


2 The model 

The Bose-Hubbard model on a lattice is defined by the Hamiltonan 

x,j/sy xsv 

where V is the set of lattice sites. The interaction U > 0 is repulsive. We restrict ourselves to nearest 
neighbour hoppings, i.e., txy = tuxy where A = {axy)x,y&v is the adjacency matrix of the lattice, i.e.. 
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ttxy = 1 if X, y are nearest neighbours, axy = 0 otherwise. We let t > 0 in order to obtain a lowest 
flat band in the systems we study. We restrict ourselves to lattices which are line graphs of planar, 
bipartite lattices. Prominent examples are the kagome lattice, the line graph of the hexagonal lattice, 
or the chequerboard lattice, the line graph of the square lattice. The usual fermionic Hubbard model 
has been studied extensively on these lattices and has ferromagnetic ground states |20| . The important 
feature of line graphs is that the spectrum of the adjacency matrix is bounded from below by —2 and 
that the eigenvalue —2 has a large degeneracy. Let |H| be the number of lattice sites and n the number 
of lattice sites per unit cell, then the degeneracy is |H|/n + 1. |H|/n is the number of states in the lowest 
band. The additional state is an element of the next band, which touches the lowest flat band. There is 
a one-to-one correspondence between the faces of the original planar, bipartite lattice and the localised 
eigenstates of A. Each face has, since the original lattice is bipartite, an even number of edges. Each edge 
is, by construction, a vertex of the line graph. A single particle state which has a constant modulus on 
these edges and an alternating sign is an eigenstate of A with eigenvalue —2. The eigenstates are linearly 
independent and thus form a basis of the eigenspace belonging to the eigenvalue —2. Fig. [^illustrates the 
procedure for the chequerboard lattice. 



square lattice 


chequerboard lattice 


localised states forming a CDW state 


Figure 1: Illustration of the construction of the chequerboard lattice as the line graph of the square lattice 
and localised states belonging to the faces. 


Let us denote the original bipartite planar lattice as G. It consists of two sub-lattices, Gi and G 2 . 
Formally, each edge of the original bipartite planar lattice can be oriented to point from Gi to G 2 . 
Further, each face can be oriented clockwise. Since the vertices x correspond to edges in G, we define for 
each face / 

1, if X belongs to the 

boundary of / and points 
into the direction of /, 


Sfx 


— 1, if X belongs to the 

boundary of / and points 


( 2 ) 


into the opposite direction 


of /, 


0, otherwise. 
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and the creation operators 
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Due to the local structure of 6 ^, it is possible to construct multi-particle ground states for the Bose- 
Hubbard Hamiltonian for densities below some critical density pc- For low densities, there exist many 
such states. Each set F of non-touch-ing faces yields a ground state \F) = fact, one can 

show that these states only form a subset of all ground states at low densities. A complete rigorous 
description of all multi-particle ground states below a certain critical density pc has been given in |19) for 
all line graphs of planar bipartite lattices. The critical density is given by a close packing of non-touching 
faces. Depending on the lattice, there may be more than one close packing. The chequerboard lattice has 
two such close packings. This procedure yields thus one or more charge density wave (CDW) states. On 
the chequerboard lattice, the critical density is pc = j. There are two such states, one of them is depicted 
in Fig. We denote the corresponding two subsets by F^ and F^- 

Above the critical density it has to be examined whether the CDW is destroyed completely or is disturbed 
merely locally. Therefore, the energies of these configurations have to be compared for a lattice filling with 
one further particle. In the case of the kagome lattice and for weak interactions, Huber et al. showed 
that the additional particles become delocalised and form a condensate. The kagome lattice has a critical 
density of pc = 5 three degenerate CDW ground states at the critical filling. As a consequence, a 
CDW state on the kagome lattice has interstitial sites, which belong to no occupied face. In this aspect, 
it differs to the chequerboard lattice, on which a CDW state has no interstitial sites. 


3 Weakly interacting model 

In this section, we study the Bose Hubbard model on the chequerboard lattice for weak interaction using a 
variational approach similar to the one used by Huber et al. | 6 ] for the kagome lattice. This section contains 
no essentially new results, we just confirm that the results in | 6 ] can be extended to the chequerboard 
lattice and presumably to any other lattice which is a line graph of a planar bipartite lattice. 

The main problem with the basis states constructed above is that these states are not orthogonal. 
To keep the property of localised states and to deal with an orthonormal basis, we follow Huber et al. 
and construct a basis of Wannier states. This allows us to project to the lowest band and to derive an 
effective Hamiltonian, which can be treated further. As a starting point, we write the hopping part of the 
Hamiltonian Q as 

Hhop = t ^ ^kk) 

k 

fzi{ki)-2 z^{ki,k 2 )\ l'bA,k\ 

\Z2{ki,k2) Zi{k2) - 2J ’ 

with 

zi {ky) = -h 2 

Z2 (fci, k2) = l + 

using /cjy = k • aj, for ly G {1,2}. The vector k belongs to the first Brillouin zone. The two vectors a,^ are 
shown in Fig. The eigenvalues are —2t and 2t [1 -|- cos(/ci) -|- cos(/c 2 )] and are plotted in Fig. [^as well. 


(4) 

(5) 

( 6 ) 
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We now construct the Wannier basis for the lower flat band. The Wannier operators are 


wj = 


hA(rj - ri)b^Aj 
j 

+w%{rj - ■ 


( 7 ) 


WAi^i) = 


(2vr)2 


1 + e 


—ik2 


[4 + 2 cos(A:i) + 2 cos(/c2)] ■ 


( 8 ) 


WB{ri) = - 


(Pk 


„*kr,: 


1 + e 


—iki 


[4 + 2cos(/i;i) +2 cos(/c 2)] 


1 5 

2 


(9) 


with the lattice vector r* = rrijuSLu and integers rriiu. The index i runs over all elementary cells of the 
lattice. The creation operators 6 ^ 4 ^^ ^ create bosons on the lattice site A or respectively B in the elementary 

cell i. This state W- involves bosons of all lattice sites, weighted with the appropriate coefficients which 
depend on the distance to the elementary cell i. These coefficients decay with increasing distance at least 
by Vlr| . 

Using the Wannier states, we write the local operators as 




■- Y. - rj)W, 

j 

+ higher band]. 


t 


( 10 ) 


Projecting onto the flat band means neglecting the contribution of the higher band. This corresponds to 
the low-energy regime for a weakly interacting system. We insert Eq. (10) into the Hamiltonian (Q and 
remove the energy offset —2tN. The largest term in the resulting Hamiltonian is an on-site repulsion. It is 
much larger than all the other terms. Following Huber et al. | 6 ], we replace it by a hard core interaction. 
The remaining terms are interactions and assisted hoppings. Due to the decay of the Wannier functions, 
the coefficients decay with increasing distance. The first leading terms are 


ffS’ = PoU ^ w'Wt 


where 




Ai^i + [^2,i + h-c.] j-Ho) 


( 11 ) 


Aia — 


h 


+ 2hwYWi 


2±ajy 


+ 2hwl 


i±ai, 




• 1/+1 


( 12 ) 
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A2,i = hW*Wi±2.. 


+ 


^idiai/dia^+1 + ^idia^q^aiy+i 


“I“^^2di2aiydia[y+i “1“ ^^2di2a^=paiy+i 


+ 


5^'''i±a. 


W^iitSai, + Wi:f2a,2 


( 13 ) 


The coefficients are 


h = 2 wa{0)^ 
h = 2wA{0fwA{-a2) 
h = 2wA{^)^WA{sii) 
h = 2wA{QfwA{aii - a2) 
h = 2wA{0fwA{-2a2) 


0.105 
-0.018 
-0.013 
0.008 
0.006 . 


(14) 


Pq projects onto states with not more than one boson in a state Wj. 

Again following Hnber et al. |6], we nse a variational ansatz for the gronnd state above pc of the form 


Z 

n n [cos(??„/ 2 ) 


0 = 1 i&Fa 


+ sin(i?o/2)w/ |0). 


(15) 


The ansatz contains two CDW states IvrOOO) = W\\^) and jOvrOO) = VbjlO) with Wa = OieFa both 
with hlling pc- Taking the expectation valne of the Hamiltonian ([TT|) yields 


- pN\^) = U 

+ Bi 

+ B 2 

+ H 3 cos(v 9 i — P 2 ) sin -di sin i ?2 


0 3 . 2 "^1 ■ 2 ^2 

2A\ sm — sm — 
2 2 


■ 2 Q • 2 ^2 ■ 2 a ■ 2 ^1 

sm tri sm — + sm 1/2 sm — 

■ 2 Q • 2 "*^1 I -20 ■ 2 ^2 

sm t/i sm-h sm i/o sm — 

2 2 


. 2 '*^1 
sm — 
2 


I • 2 "*^2 

+ sm — 





2 


• 2 ^1 I ■ 2 "*^2 
sm-h sm — 


The coefficients are combinations of the interaction coefficients /j. They have the following valnes 


Al « 0.0527 
Hi « -0.0138 


B 2 ^ -0.0002 
H 3 PS -0.0072. 


(16) 


(17) 


Eq. (16) is minimised nnmerically by varying the parameters i?i, 2 ) Pi — P 2 and /r; the latter in order to 
keep the particle nnmber hxed. The density p is varied between the critical hlling factor Pc= \ and p = \- 
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Starting point is the CDW state W"||0) at pc- This enables to determine quantitatively the expectation 
value of and by analogy to jB] the order parameter of the CDW 


and of the superfluid 


V’CDW = -—^ +-+-^—- 

{^wlWi + WlW2\^) 

cos — cos 'd 2 
cos + cos — 2 ’ 


V^sf 


2 

N 




1 

2 


sin '!?2 


(18) 


(19) 


The results are plotted as a function of the density in Fig. The order parameter of the CDW 
decreases continuously from 1.0 by increasing the density above the critical filling factor, and conversely 
the superfluid order parameter grows from 0.001 on. For densities higher than p ~ 0.36 the order parameter 
of the CDW abruptly vanishes and the superfluid phase starts to decay. 

As the set of variational states is restricted to hard-core bosons we exclude double occupancy. Hence, 
for large densities this approach is not a sufficient choice for calculations of ground states and will not 
provide reliable results. Furthermore, the truncation of the Wannier state - regarding only leading terms 
- also requires adequate low lattice fillings. Huber et al. observed a similar behaviour which might 
also confirm that we reach a limit of our model at this density regime and do not observe any physical 
response. 

We tested the stability of the approximation in Eq. 0 by taking all coefficients into account up to 
O.OOID. We observe a similar result. The onset of the decline for the superfluid order parameter remains 
more or less unchanged with 0.37. However, the density for which the order parameter of the CDW will 
be zero is shifted to 0.47. In general, this indicates, nevertheless, the stability of the truncation procedure. 

To conclude, the results for the chequerboard lattice are similar to the results for the kagome lattice 
derived by Huber et al. j6], although, as pointed out above, there are no interstitial sites on the chequer¬ 
board lattice. The essential point is, indeed, that the additional particles are delocalised. These particles 
have a high mobility and this is energetically favourable to localised states, as pointed out by Huber et al. 


4 Strongly interacting model 

If U is much larger than t, states where a site is occupied with more than one boson are energetically 
unfavourable. We take the limit U —)• oo, i.e. hard-core bosons. For densities p < 1 this means that we 
project onto the subspace of multi-particle states where each site is occupied with at most one boson. Let 
P <1 be the projector onto this subspace. Then, the Hamiltonian reads 

H = tP<i E (20) 

x,y£V 


4.1 The chequerboard lattice 

The CDW is the exact ground state at the critical density jl9| . Therefore, as in the case of small U, we 
use the CDW at the critical density and add one further boson. A first variational ansatz for the ground 
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state is thus 


\cP) = P<iY,^A\F,^) 

xev 

= Y1 -jrfi Y1 ^xSfc^'blbl- 

f^po V l/l x,x'ef,x^x' 

• n ( 21 ) 

/'6Ff\{/} 

The variation was done on a 20x20 cutout of the chequerboard lattice. For reasons of numerical stability, 
we choose a symmetric sub-unit of 16 lattice sites containing five faces on which (j)x are varied freely and 
we let (fix = (fp outside that region. The resulting minimal energy is —{2\Ff \ + 0.67)f and belongs to 
a state where the additional particle is localised on one of the faces of the sub-unit which are already 
occupied. It seems energetically favourable to disturb only one particle instead of four when putting the 
additional particle on a not occupied face. It seems also energetically favourable to have the additional 
particle localised and not distributed over many faces. 

Since the additional particle is localised on an already occupied face, a second variational ansatz for the 
ground state is 

\u,f)= Y. Uee'hlhl,bf\Ff), (22) 

e,e'e/eFf 

for an arbitrary face of Tf. Since all faces are independent, this variational problem is equivalent to solving 
the problem for two hard-core bosons on a square. The resulting energy is —2^/2t for the two hard-core 
bosons and therefore —{2\Ff\ — 2 + 2^/2)t = —{2\Ff \ -|- 0.83)f. If this was the exact ground state, it would 
be degenerate with a degeneracy As a consequence of this variational ansatz, we obtain an exact 

upper limit for the ground state energy per particle eo(/9) for hard-core bosons on the chequerboard lattice 
at density p to be 

eo(/o) <-[2pc + (2\/2 - 2)(p - Pc)]T (23) 

which holds for Pc ^ P < ‘^Pc- The number of states at or below that energy must be equal or above the 
number of occupied faces. 

As a further test of the ansatz, we look at the next non-trivial sub-unit, depicted in Fig. The exact 
ground state with four particles is the state where the four faces in the corners are occupied, depicted as 
circles in Fig. We now put a fifth boson into the sub-unit. Following the argument above, we conclude 
that there are at least four ground states with an energy below — (6 -|- 2^/2)t = —8.83t. A numerical 
diagonalisation of the system yields, indeed, exactly four states with energies below that value. These 
energies are —9.016t, —9.022t (two-fold degenerate), and —9.081t. The corresponding eigenstates have a 
large overlap with combinations of the four variational states. They thus show a similar behaviour as it 
was observed in different one-dimensional models dsnniiiB]. Pairs of bosons are formed. 

Since this sub-unit shares no sites with other occupied faces, we can use this result as a new upper 
variational limit for densities between pc and |pc) 

eo{p) < -[2pc +1.081 {p - pc)]t. (24) 


4.2 The kagome lattice 

Whereas at weak interactions, physics of hard-core bosons on the chequerboard and the kagome lattice 
are similar, there may be a difference at strong interactions. On the kagome lattice a CDW ground state 



at the critical density has interstitial sites, whereas on the chequerboard lattice there are no interstitial 
sites. 

On the kagome lattice, the smallest building block is a hexagon. One particle on the hexagon has an 
energy of —2t. If we put two hard-core bosons on a hexagon, the energy is —2^/3t. 

The smallest non-trivial sub-unit of the kagome lattice consists of one hexagon, surrounded by six 
triangles and three non-touching hexagons, as depicted in Fig. 

Three hard-core bosons on this sub-unit are placed on the three outer hexagons; the corresponding 
state is the only ground state on this sub-unit with energy —6t. This state has three interstitial sites. The 
ground state with four particles on this sub-unit must have an energy below —(4 -|- 2^/3)t = —7A6t and 
following our variational argument, there are at least three eigenstates below this energy. 

A numerical diagonalisation of the sub-unit with four particles shows that there are three almost de¬ 
generate states with energies between —7.55t and —7.561 and one state, the ground state, with an energy 
of —7.66t. The ground state has a large expectation value of the particle number on the three interstitial 
sites of the inner hexagon, whereas this expectation value is small on the other three low eigenstates. The 
overlap of the ground state, with one of the states where two particles are on a hexagon, is only 0.085. The 
other three exact eigenstates with energies below —{4:+2^/3)t have a high overlap with linear combinations 
of the three states with two particles on a hexagon. The overlap lies above 0.95. This shows, that linear 
combinations of the variational states, where the additional hard-core boson is put on one of the occupied 
hexagons, describe these three low energy states well. However, in the true ground state of the sub-unit 
the additional particle moves on the interstitial sites and is not localised on one of the occupied hexagons. 

This argument shows that interstitial sites on the kagome lattice are important if hard-core bosons are 
added to one of the ground states at the critical Filing. This suggests that the argument by Huber et al. 

may remain correct in the regime of strong interaction. The additional bosons move mainly on the 
interstitial sites and form a condensate at densities slightly above the critical density. 

4.3 Other line graphs 

The same arguments apply to line graphs of other bipartite plane lattices. If the faces of the underlying 
bipartite plane lattice form a two-colour map, the CDW states are at most two-fold degenerate. Since 
each edge belongs to two faces, there are no interstitial sites. On the other hand, if three or more colours 
are needed to colour the faces of the underlying bipartite plane lattice, the CDW states have interstitial 
sites. This is true for inhnite lattices and also for finite sub-units. In the latter case the outer faces has to 
be coloured as well. Therefore, the sub-units should be chosen such that they need the same number of 
colours as the original lattice. 

An example for the case with more than two colours is the so called truncated square tiling or 4.8^- 
tiling, a bipartite plane lattice formed of octagons and squares. Hard-core bosons on its line graph have a 
unique CDW ground state with one boson on each square. There are interstitial sites stemming from the 
additional edges of the octagons. The smallest non-trivial sub-unit consists of one octagon surrounded by 
four non-touching squares. Putting hve hard-core bosons on this sub-unit, we expect at least four ground 
states with an energy less than —{6 + 2^/2)t, corresponding to the four states where the additional particle 
is put on one of the already occupied squares. A numerical diagonalisation of the sub-unit yields eight 
eigenstates with an energies less than —{6 + 2\/2)t. As for the kagome lattice, the ground state has a large 
expectation value of the particle number on the four interstitial sites. Let us mention that in this case it 
might even be favourable to put two bosons on the interstitial site, because two hard-core bosons on the 
octagon have an energy of —\/l0 -|- 2\/5t = —3.804t. This state contributes to the low lying eigenstates 
as well. The analysis on this sub-unit thus shows that the ground state is dominated by interstitial sites. 

An example where the faces of the underlying bipartite plane lattice form a two-colour map, like for the 
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square lattice, can be constructed from the triangular lattice. We put an additional lattice site on each 
edge of the triangular lattice. The resulting lattice is a bipartite plane lattice. Its line graph has hexagons 
connected via complete graphs with six vertices, similar to the kagome lattice, where the hexagons are 
connected via complete graphs with three vertices. Two particles on a hexagon have the energy —2\/3. 
Let us take a sub-unit with four hexagons, one in the middle surrounded by three non-touching hexagons. 
Since there are three non-touching hexagons, we expect at least three eigenstates with energies below 
— (4 -|- 2y/i)t. The exact diagonalisation yields four such eigenstates. The three lowest states have, again, 
a large overlap between 0.66 and 0.88 with linear combinations of the states where the additional particle is 
put on a hexagon. The fourth eigenstate appears because of the large boundary of the sub-unit. Indeed, if 
we put two particles on the inner hexagon and two particles on the ring of 12 sites forming the boundary, 
we get a variational state with an energy of —7.36f, only slightly above —(4 -|- 2^/3)t. This additional 
variational state has thus two bound pairs of bosons, in contrast to the three others with one bound pair, 
and a large overlap with two of the four low lying eigenstates. Thus, in the ground state and in the three 
low lying eigenstates of this sub-unit, bound pairs of bosons are formed as for the chequerboard lattice. 

There exist one-dimensional systems where these arguments apply as well. Consider for instance the 
line graph of the square chain, a lattice of two chains of corner sharing triangles. The square chain needs 
indeed three colours because the two outer faces cannot be coloured with the same colour as the squares. 
Thus, the CDW states have interstitial sites. A numerical diagonalisation shows indeed that for hard¬ 
core bosons, additional particles move on these sites. Let us mention that this one-dimensional lattice is 
sometimes called kagome-chain, see e.g. |T^ 

5 Conclusion and outlook 

We showed that weakly interacting bosons on the chequerboard lattice behave similarly to weakly inter¬ 
acting bosons on the kagome lattice. We essentially confirmed the results obtained by Huber et al. for the 
kagome lattice. Above the critical density the additional particles are delocalised and form a condensate. 

For hard-core bosons, the situation is different: our results indicate that on the chequerboard lattice 
the additional particles form pairs of two bound bosons. This is similar to the findings of Phillips et al. 
for the saw-tooth chain )18j . In contrast to that, on the kagome lattice the additional particles move on 
interstitial sites as for weak interaction. 

The kagome lattice and the chequerboard lattice are just examples of two larger classes of two-dimensional 
lattices with flat bands. The kagome lattice is a prototype of line graphs of bipartite plane lattices where 
the faces cannot be coloured with two colours. Instead, three or four colours are needed. For all lattices 
in this class, the CDW states have interstitial sites and we expect that additional particles will move on 
these sites, even for strong interactions. The chequerboard lattice is a prototype of line graphs of bipartite 
plane lattices where the faces form a two-colour map. In that case, since each edge of the original bipartite 
plane lattice belongs to exactly two faces, the CDW states have no interstitial sites and pairs of bound 
particles are formed for strong interactions. 

There are several questions which remain open. For instance, one may easily construct line graphs 
of three dimensional lattices like the line graph of the diamond lattice, an analog to the kagome lattice 
in three dimensions, or the line graph of the simple cubic lattice. Partly, such lattices can be found in 
nature. The line graph of the diamond lattice is the octahedral sublattice of a spinel. Being line graphs 
of bipartite graphs, these lattices have a lowest flat band and localised single particle eigenstates in that 
band. As in two dimensions, it is possible to construct strictly localised states. As in two dimensions, 
they are localised on the elementary cycles, squares or hexagons in the two examples. But unfortunately, 
a complete description of all multi-particle ground states at or below the critical density is not known. 
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The arguments in |19| do not apply here. Therefore, a discussion of these lattices above the critical hlling 
is out of reach today. 

Let us mention that hard-core bosons on graphs is a topic that has recently been discussed in (algebraic) 
graph theory in connection with the isomorphism problem j22l [23l [2^ |25] . Spectral properties of hard¬ 
core bosons on a graph can reveal certain properties of the graph which are otherwise difficult to obtain. 
There is no direct relation between the current paper and these investigations, which mainly deal with 
strongly regular graphs. Nevertheless, it is interesting to see that here as well a certain spectral property 
of hard-core bosons is connected to a graph theoretical problem, here the colouring of the faces. 
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Figure 2: The band structure of the chequerboard lattice in the first Brillouin zone (right) and the choice 
of the unit cell (left) with the two lattice vectors ai ^2 and the two lattice types A and B. 



Figure 3: The expectation value of [cf. Eq. (16)] is plotted on the left y-axis for the parameters 

obtained by solving the variational problem. Eq. (18) and (19) are likewise depicted on the 
right y-axis for the order parameter of the CDW and the superfluid, respectively. 
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Figure 4: Sub-unit of the chequerboard lattice. The cycles denote the occupied faces. In the exact ground 
state with four particles, each face is occupied by one boson. In the variational states (22) with 
five bosons one of the faces is occupied by two bosons. 



Figure 5: Sub-unit of the kagome lattice. The cycles denote the occupied 
with three particles, each face is occupied by one boson. In the 
bosons one of the faces is occupied by two bosons. 


faces. In the exact ground state 
variational states (22) with four 
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